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Context. Various observations of magnetic fields in the intra-cluster medium (ICM), most of the time restricted to cluster cores, point towards 
a field strength of a few /jG (synchrotron radiation from radio relics and radio halos, inverse Compton radiation in X-rays and Faraday rotation 
measure of polarised background sources). Both the origin and the spatial structure of galaxy cluster magnetic fields are still under debate. In 
particular, the radial profile of the magnetic field, from the core of clusters to their outskirts, is important for cosmic ray propagation within the 
cosmic web. 

Aims. In this letter, we highlight the importance of cooling processes in amplifying the magnetic field in the core of galaxy clusters up to one 
order of magnitude above the typical amplification obtained for a purely adiabatic evolution. 

Methods. We performed a "zoom" cosmological simulation of a 3 keV cluster, including dark matter and gas dynamics, atomic cooling, UV 
heating, and star formation using the newly developed MHD solver in the AMR code RAMSES. 

Results. Magnetic field amplification proceeds mainly through gravitational contraction. Shearing motions due to turbulence provide additional 
amplification in the outskirts of the cluster, while magnetic reconnection during mergers causes magnetic field dissipation in the core. 
Conclusions. Cooling processes have a strong impact on the magnetic field structure in the cluster. First, due to the sharp rise of the gas density 
in the centre, gravitational amplification is significantly higher, when compared to the non-radiative run. Second, cooling processes cause 
shearing motions to be much stronger in the core than in the adiabatic case, leading to additional field amplification and no significant magnetic 
reconnection. Cooling processes are therefore essential for determining the magnetic field profile in galaxy clusters. 
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1. Introduction 

Cl usters of galaxies are kn own to be magnetised (see review 
by iGovoni & Ferettil I2004I) . The existence of magnetic fields 
has been determined either by direct methods like diffuse syn- 
chrotron radio or inverse Compton hard X-ray emission or by 
indirect methods like Faraday rotation measures (RM). They all 
suggest that ytG fields lie in central regions of galaxy clusters 
up to several 10 fiG magnetic fields in large cooling flows clus- 
ters. Magnetic field strength can differ from one method to the 
next Direct methods usually capture large-scale fields averaged 
over large volumes, while RM are derived from the analysis 
of background point sources and are thus sensitive to small- 
scale variations (cold filaments, shear flows, shocks, galaxy 
stripping, galaxy winds, etc.). To shed light on the magnetic 
topology found in cosmic structures, it is of great interest to 
perform direct, self-consistent numerical simulations of galaxy 
clusters. Magnetic fields in clusters are also important for de- 
termining the deflection angle of ultra-high energy cosmic rays, 
since they probably host the source of these cosmic rays (s ee 
LemoineM2005 : Kot era & Lemoinel2008l : Glob us et al 120071) . 
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Simulations of galaxy clusters with magnetic fields have 
been p erformed using smooth par t icle h y drodynamics (SPH ) 



codes dDolag et all 1 19991 iDola 
grid-based cod e s (iRoettiger et 
Sigl et all l2004t lAsai et all 120071) 



ag 
al. 



20001; iDolag et alJ |2005l) 
19991: iMiniati et al.l 12001 



a nd ada ptive mesh refine- 



ment (AMR) codes (iBruggen et al. 2005b. using both cos 
mological simu l ations dDolag et alJ 1999 :" lM~iniati et ak 2001 



Sigl et al.1 120041: IDolag et al.l 120051) an d idealised simulations 
(IRoettiger et alJI 19991: lAsai et al.ll2007l) . In this letter, we report 



the first cosmological simulation with AMR to include atomic 
cooling, UV heating, and star formation physics, with a full 
treatment of the ideal MHD equations. We have also performed 
a reference adiabatic run to compare our results with previous 
works and to point out the differences with the radiative case. 

2. Simulations 

We performed a "zoom" cosmological simul ation of a galax y 



cluster using the AMR code RAMSES dTeyssierl l2002h . 



Gas dynamics was computed using a second -order unsplit 



Godunov scheme for the ide al MHD equations dTeyssier et al 



20061: iFromang et alJ 120061) . while collisionless dark matter 
particles are evolved using a particle-mesh solver. Gas cool- 
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ing and heating were taken into account as source terms 
in the energy equation. The cooling and heating functions 
were computed for a pr imordial H and He plasma, using the 
Haardt & Madaul (y.996) background model. Radiative losses 
lead to the formation of high density, (low temperature) re- 
gions, where stars are allowed to form according to a Schmidt 
law: p* = ep/tff if p > po- The density threshold for star forma- 
tion was set to po = 10 5 Q£,p c (z). The star formation efficiency 
was set to e = 5%. The simulation comoving box length was 
chosen equal to 80 tr 1 Mpc with a ACDM cosmology with 
Q m = 0.3, Q A = 0.7, £l b = 0.045, H = 70km.s~ 1 .Mpc- 1 , and 
cr 8 = 0.9. A spherical region of radius 12.5 h -1 Mpc around 
our simulated cluster was defined as our high-resolution region, 
with an effective resolution of 512 3 . A coarser grid with an ef- 
fective resolution of 256 3 was used to cover the inner 40 hr 1 
Mpc, and finally an even coarser 128 3 grid was used to cover 
the whole box. The mass of dark matter particles on each coarse 
grid are 2.9 x 1O IO M , 3.6 x 1O 9 M , and 4.5 x 1O 8 M . Only 
the finest grid was allowed to trigger new refinements during 
the course of the simulation, up to 7 additional levels of AMR 
cells leading to a maximum resolution of 1.2h _1 kpc. We used 
a quasi-Lagrangian criterion: each cell is individually refined 
if the number of dark matter particles exceeds 8, or if the bary- 
onic mass exceeds 8 times the initial high-resolution mass res- 
olution. We solved the full set of ideal MHD equations using 
a new scheme based on a Goduno v implementat i on of con- 
strained transport and presented in iTevssier et al. I d2006l) and 
Fromang et al. J 20061) and we u sed the HLLD Riemann solver 
from lMiyoshi & Kusanol (120051) . The comoving magnetic field 



was set initially to a const ant value, B z ^ 10 11 G, as sug- 



gested in lDolag et alj d2005l) . to reproduce the pG fields in clus- 
ter cores. With these parameters in the course of the simulation, 
the plasma B = P ga s/Pmag never decreased below 1000: the dy- 
namical effect of the magnetic field can therefore be considered 
as negligible, even in the core of our simulated cluster. 

3. Results 

Figure Q] shows the column density distribution of the gas and 
the mean magnetic field for the cooling run at z = 0. Magnetic 
field amplitudes are strongly correlated with the density dis- 
tribution with mass-averaged values of B ~ 10~'/iG in the 
cluster core, galaxies with B ~ 10~ 2 pG in satellite clumps, 
and B ~ 10 ~ 3 pG in filaments. We define the virial mass as 
M200 = 200 x 47r/3p e 7?2 00 , where p c is the critical density. 
For the adiabatic simulation, we found the following prop- 



erties for our cluster at z = 0: R'^ 



= 1 h-'Mpc, M%h 



2.7 x 10 14 h _I M o , and Tf - 3.4 keV. In the radiative case, 
we obtained /?200 - 1 1 h~'Mpc, M 2 qq - 3.5 x 10 14 h _1 M o , and 
Tx - 5.1 keV. The magnetic field amplification of a collaps- 
ing three-dimensional gas sphere with infinite conductivity is 
given by B oc p 2 ^ 3 , for the magnetic flux to be conserved. Thus 
a subsequent increase (or decrease) in the magnetic field with 
respect to this purely gravitational amplification should reveal 
other amplification or dissipation mechanisms. Figure [2] shows 
the mass-weighted histogram of the radiative simulation in the 
p-\B\ plane, within 2 virial radii around the cluster. For densities 
lower than 10 4 p, the magnetic field amplification is one order of 



Fig. 1. Logarithm of the column density map of the gas in 
units of the mean baryons density (upper panel), and of the 
mass-averaged magnetic field amplitude in units of pG (bot- 
tom panel) at z — 0. 



magnitude h igher than for pure gravitational compression. As 
discussed in iDolag et al.l (12005b . this is likely due to shearing 
motions in the cluster atmosphere, caused by turbulence and 
frequent mergers. At higher densities, the radiative run diverges 
strongly from the adiabatic case. The gravitational compression 
due to the cooling flow provides additional field amplification 
in the high-density tail. As seen below, cooling also provides 
a sustained turbulent regime in the core and the corresponding 
additional field amplification. Bas ed on the Zel'dovich a pproxi- 
mation of gravitational dynamics. lKing & Coles (2006) predict 
that a cosmological magnetic field should evolve as a B oc p 87 , 
due to anisotropic collapse in a Gaussian random field. This 
compares favourably with the low-density part of our simula- 
tion (see Fig. [2]). 

At higher density, the situation is more complex in the 
core of our simulated cluster. In the adiabatic case, the mag- 
netic field amplitude decreases below the expected value for 
pure compression. The mean magnetic field is of the order of 
10~ 2 /zG, far below the typical values of magnetic field ampli- 
tude in observed cluster cores. It is also apparent in the mag- 
netic field profile plotted in Fig. [3] for which a dip in the field 
strength is visible in the cluster core. While shear flows are able 
to sustain additional magnetic amplification in the outer parts 
of the cluster r > 150 fr'kpc, we find the opposite behaviour in 
the cluster core where magnetic field is dissipated. According 
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Fig. 2. Histogram of the mass fraction for the cooling run as a 
function of the normalised density and the magnetic amplitude 
at z — 0. The black solid line is the p 2 ^ 3 collapse amplification, 
the averaged magnetic field as a function of the normalised den- 
sity for the cooling case (dashed line) and the adiabatic case 
(dotted line) are also plotted. 




Fig. 3. Mean magnetic field as a function of the normalised ra- 
dius. The black lines are the p 2/i collapse amplification, the 
green lines denote the run values for the adiabatic (dotted) and 
cooling (solid) runs at z — 0. 



to iRoettiger et al.l d 19991) . this effect stems from magnetic re- 
connection occurring during merger events. Since we are not 
considering any microscopic process here and since the MHD 
turbulence scale is far below our numerical resolution, this dis- 
sipation is caused by our numerical scheme that captures the 
weak solution of the ideal MHD equations. Although mag- 
netic reconnection probably occurs in nature within converg- 
ing flows, the exact amplitude of the dissipation is likely to de- 
pend on the microphysics. In the present numerical approach, 
results should depend strongly on the spatial resolution and on 
the numerical scheme used. It is however interesting to anal- 
yse the effect of cooling in this respect. As can be seen in 
Fig. [3j magnetic dissipation in the cluster core is suppressed 
and magnetic amplification now proceeds in the same way as 
in the outer parts, with gravitational compression and shearing 
motions. Only in the very centre (below 3 fr'kpc, close to the 
resolution limit) do we see magnetic dissipation again. 




Fig. 4. Mean radial velocity (black), radial velocity dispersion 
(blue), and sound speed (red) for the adiabatic (solid) and the 
cooling (dashed) runs at z — in units of V200 for the adiabatic 
case. 

To illustrate this point further, we show in Fig. [4] velocity 
profiles in the adiabatic and in the cooling cases. The radial ve- 
locity dispersion is a signature of turbulent motions, so it is no 
surprise to see that a strong velocity dispersion at a given ra- 
dius corresponds to an excess of field amplification at the same 
radius (see Fig. [3). In the adiabatic case, turbulence is dissi- 
pated in the core, and magnetic dissipation occurs. In the cool- 
ing case, gravitational contraction resumes, as well as shearing 
motions, so that magnetic amplification is now more active in 
the core. The RM maps at z = of the adiabatic simulation and 
the cooling simulation shown in Fig. [5] differ strongly. In the 
adiabatic case, we obtain RM values of 30 rad.irT 2 in the clus- 
ter centre, whereas in the radiative case, we reach 1000 rad.irT 2 
in the very centre (r < 10 fr'kpc) and obtain RM values of 
lOOrad.rrT 2 in the core. The RM results of the cooling sim- 
ulation are marginally consistent with the IClarke et al.l (1200 ll) 
sample (200 rad.irT 2 ) a nd also with t he ma ximum values found 
in the cluster sample of lTavlor et al.l (120021) (up to 1800 rad.irT 2 



for the hot gas cluster). The magneti c field in the outer p art of 
our cluster is on the weak side of the Clarke et al. ( 2001b sam- 



ple. It could result from differences in cluster mass or from our 
somewhat arbitrary choice of initial magnetic field. Note that 
we can increase the initial amplitude up to 10 times while re- 
maining in the linear regime (J3 > 10), so that all our results can 
be re-normalised by this factor. 

4. Conclusion and discussion 

There are noticeable differences in the magnetic field character- 
istics between a galaxy cluster with an adiabatic evolution and 
a galaxy cluster with radiative cooling: the average magnetic 
field in the cluster core is significantly higher when a cool- 
ing flow is present, due to additional gravitational compression 
but also an increased level of turbulence in the core-driving 
shearing motions. The main consequence is that Faraday RM 
simulated maps agree more with observations in the cooling 
case, if the initial comoving magnetic field value is taken to 
be 10" G. In low-density regions, however, the magnetic 
field evolution in the radiative run is very close to the adia- 
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Fig. 5. |RM| maps colour scale of the cluster core for the adia- 
batic run (left panel) and the cooling run (right panel) at z = 0. 



batic case since dynamical properties of the gas turbulence are 
nearly the same. Note that the turbulent cascade is far from 
being fully resolved in our simulation, especially outside the 
cluster core where our resolution barely reaches lOIr'kpc, a 
muc h larger scale than the dissipation scale o f MHD turbulence 
(e.g. Brunetti & Lazarianl200 7t ljonesll2007h . For these reasons 
our results should only be considered as lower estimates for the 

magnetic field strength. 

As was already discussed in iRoetti geret al.ldl999h in the 
context of adiabatic simulations of idealised mergers, magnetic 
reconnection might be responsible for field dissipation in the 
cluster core. Since we are dealing with ideal MHD, magnetic 
dissipation occurs at the numerical level, so that we should 
be affected to some extent by the effect of numerical resolu- 
tion. Moreover, this underlines the importance of the choice 
of the numerical code used, especially when one considers 
the fundamental differences between grid-based a nd particle- 
based codes. Using an MHD version of GADGET. iDolag et al 
(120051) find a much larger magnetic field strength in their adi- 
abatic run for a subset of their simulated particles. They re- 
port mean field values one order of magnitude higher than 
median field values, apparently contradicting our present re- 
sult. Magnetic dissipation therefore appears much less efficient 
in the SPH case. The Lagrangian nature of this technique is 
radically different from grid-based codes. Since both codes 
lead to different magnetic diffusion terms, a proper compari- 
son between numerical techniques on this stake would be very 



helpful to the commu nity. Using the grid-based code ZEUS, 
Roettiger et al] (119991) report very similar results to ours, with 
strong field dissipation occurring in the converging part of the 
flow. One interesting outcome of the present work is that radia- 
tive cooling drastically changes the effect of magnetic recon- 
nection, since turbulence and gravitational compression easily 
counterbalance the associated dissipation. 

Nevertheless, this large additional amount of magnetic en- 
ergy in the core of cooling clusters is a crucial step in deter- 
mining the structure of the cosmological magnetic field. It has 
a direct consequence on the propagation of high-energy cos- 
mic rays in the universe. Since we have no direct observations 
of magnetic fields outside of cluster cores, only cosmological 
numerical simulations can address this problem. Their weak- 
ness is that the initial magnetic field value must be normalised 
a posteriori in order to fit the observed values. We have shown 
that, for the same seed field, the final magnetic field strength in 
a cooling cluster core is one order of magnitude more than in 
the adiabatic case. Cooling processes are therefore important if 
one wants to describe the proper evolution of magnetic fields 
in the Universe. 
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